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Abstract 

Electric-field-driven transport of electronic fluids in metallic glasses as well as three- 
dimensional amorphous superconductors are investigated by using the verified approach 
which has been successfully adopted to study the critical transport of glassy solid helium 
in very low temperature environment. The critical temperatures related to the nearly fric- 
tionless transport of electronic fluids were found to be directly relevant to the superconduct- 
ing temperature of amorphous superconductors after selecting specific activation energies. 
Our results imply that optimal shear-thinning is an effective way to reach high-temperature 
charged superfluidity or superconductivity. 
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1 Introduction 

Analogies between superconductivity (charged superfluidity) and superfluidity in liquid [1-2] and 
solid helium (either supersolid [3] or superglass [4-6]) has been noted by some researchers. The 
particles (charged electrons as well as atoms) in a superfluid are in the same quantum state, 
so they move coherently and result in dissipationless mass flow. Quite recently a glassy state 
of supersolid helium (solid 4 He : ability to flow without resistance) has been reported [5]. We 
have made some contributions to the possible transport of glassy solid 4 He in confined domain 
under rather low temperature environment (cf. Ref. 6 for the comparison with available relevant 
experiments). Our present interest is to borrow the above mentioned analogy to investigate the 
nearly frictionless states for transport of the electronic liquid (presumed to be glassy or amor- 
phous) in amorphous superconductors [7-10] which is closely related to charged frictionless or 
superfluidity states. The approach adopted in [6] (considering transport of glassy fluids) will be 
extended by varying the forcing to be an electric field in this paper. 

Note that one prominent difference between the fluid transport in microdomain and those in 
macrodomain is the strong fluid-wall interactions observed in microconduits. For example, as 
the microconduit size decreases, the surface-to-volume ratio increases. Therefore, various prop- 
erties of the walls, such as surface roughness, greatly affect the fluid motions in microconduits. 
In this paper, we adopt the verified Eyring model [6,11-13] to study the nearly frictionless trans- 
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port of (glassy) electronic fluids within corrugated microannuli. Eyring model treats transport 
of shear-thinning (which means once material is subjected to high shear rates then the viscosity 
diminishes with increasing shear rate) matter quite successfully [6,11-13]. To obtain the law 
of shear-thinning fluids for explaining the too rapid annealing at the earliest time, because the 
relaxation at the beginning was steeper than could be explained by the bimolecular law, a hy- 
perbolic sine law between the shear (strain) rate : 7 and (large) shear stress : r was proposed 
and the close agreement with experimental data was obtained [11-13]. This model has sound 
physical foundation from the thermal activation process (Eyring [12] already considered a kind 
of (quantum) tunneling which relates to the matter rearranging by surmounting a potential 
energy barrier). With this model we can associate the (shear-thinning) electronic fluid with 
the momentum transfer between neighboring clusters on the microscopic scale and reveals the 
detailed interaction in the relaxation of flow with dissipation (the momentum transfer depends 
on the activation shear volume, which is associated with the center distance between (charged) 
particles and is proportional to UbT/tq (T is temperature in Kelvin, and To a constant with the 
dimension of stress). Thus, this model could be applied to study transport of (shear-thinning as 
well as amorphous or glassy) electronic fluids in microdomain which is of particular interest to 
researchers in condensed matter physics (considering amorphous superconductors [7-10]). We 
select the microannular geometry as it is more relevant to the experimental environment con- 
sidering the measuring conductivity of samples of amorphous superconductors [7-10] (transport 
of shell- like electronic fluids). 

To consider the more realistic but complicated boundary conditions in the walls of microannu- 
lus, however, we will use the boundary perturbation technique [13-14] to handle the presumed 
wavy-roughness along the walls of microannuli. To obtain the analytical and approximate solu- 
tions, here, the roughness is only introduced in the radial or transverse direction. The relevant 
boundary conditions along the wavy-rough surfaces will be prescribed below. We shall describe 
our approach after this section : Introduction with the focus upon the boundary perturbation 
method. The expression of the electric-field-driven volume flow rate of the glassy (or amorphous) 
electronic fluid is then demonstrated at the end. We noticed that the conductivity measurements 
reported in [9] and [10] are primary in Mg-based and Fe^Mi^-based metallic glasses. Finally, 
we will illustrate our results into four figures and give discussions therein. Our results show 
that some specific activation energies and activation volumes are crucial to the nearly friction- 
less states of the transport of the electronic fluids at rather low temperature regime. The good 
qualitative as well as quantitative (e.g., critical temperature or T c ) comparison with experiments 
confirms our approaches. Our results imply that optimal shear-thinning is an effective way to 
reach high-temperature charged superfluidity or superconductivity. 
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2 Formulations 

Researchers have been interested in the question of how amorphous material responds to an 
external mechanical load. External loads cause (electronic) liquids to flow, in Newtonian or 
various types of non-Newtonian flows. Glassy materials, composed of polymers, metals, or 
ceramics, can deform under mechanical loads, and the nature of the response to loads often 
dictates the choice of material in various industrial applications. The nature of all of these 
responses depends on both the temperature and loading rate. As described by Eyring [12], 
mechanical loading lowers energy barriers, thus facilitating progress over the barrier by random 
thermal fluctuations. The Eyring model approximates the loading dependence of the barrier 
height as linear. The Eyring model, with this linear barrier height dependence on load, has been 
used over a large fraction of the last century to describe the response of a wide range of systems 
and underlies modern approaches to sheared glasses [11-12]. To the best knowledge of the author, 
the simplest model that makes a prediction for the rate and temperature dependence of shear 
yielding is the rate-state Eyring model of stress-biased thermal activation [11-12]. Structural 
rearrangement is associated with a single energy barrier E that is lowered or raised linearly by 
an applied stress a. 

R± = voexp[-E/(k B T)]exp[±(TV*/(kBT)], 

where ks is the Boltzmann constant, vq is an attempt frequency and V* is a constant called the 
activation volume. In glasses, the transition rates are negligible at zero stress. Thus, at finite 
stress one needs to consider only the rate R+ of transitions in the direction aided by stress. 
The linear dependence will always correctly describe small changes in the barrier height, since it 
is simply the first term in the Taylor expansion of the barrier height as a function of load. It is 
thus appropriate when the barrier height changes only slightly before the system escapes the local 
energy minimum. This situation occurs at higher temperatures; for example, Newtonian flow is 
obtained in the Eyring model in the limit where the system experiences only small changes in the 
barrier height before thermally escaping the energy minimum. As the temperature decreases, 
larger changes in the barrier height occur before the system escapes the energy minimum (giving 
rise to, for example, non-Newtonian flow). In this regime, the linear dependence is not necessarily 
appropriate, and can lead to inaccurate modeling. To be precise, at low shear rates (7 < j c ), 
the system behaves as a power law shear-thinning material while, at high shear rates, the stress 
varies afhnely with the shear rate. These two regimes correspond to two stable branches of 
stationary states, for which data obtained by imposing either a or 7 exactly superpose. The 
transition from the lower branch to the higher branch occurs through a stable hysteretic loop 
in a stress-controlled experiment. 

The motion (of composite particles) is represented in the configuration space; on the potential 
surface the stable particles are in the valleys, which are connected by a pass that leads through 
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the saddle point. A (composite) particle at the saddle point is in the transition (activated) 
state. Under the action of an applied stress the forward velocity of a flow unit is the net number 
of times it moves forward, multiplied by the distance it jumps. Eyring proposed a specific 
microscopic model of the amorphous structure and a mechanism of deformation kinetics [11-12]. 
With reference to this idea, this mechanism results in a (shear) strain rate given by 
. V h k B T —AE . V h r 

where 

Vh = A2A3A, V m = A2A3A1, 

Ai is the perpendicular distance between two neighboring layers of (composite) particles sliding 
past each other, A is the average distance between equilibrium positions in the direction of 
motion, A2 is the distance between neighboring (composite) particles in this same direction 
(which may or may not equal A), A3 is the (composite) particle to (composite) particle distance 
in the plane normal to the direction of motion, and r is the local applied stress, AE is the 
activation energy, h is the Planck constant, fee is the Boltzmann constant, T is the temperature, 
Vh is the activation volume for the microscopic event [12]. The deformation kinetics of the chain 
(composite particles) is envisaged as the propagation of kinks in the particles into available 
holes. In order for the motion of the kink to result in a flow, it must be raised (energised) into 
the activated state and pass over the saddle point. This was the earliest microscopic theory of 
yield behaviour in amorphous materials, and Eyring presented a theoretical framework which 
formed the basis of many subsequent considerations. 
Solving Eqn. (1) for r, one obtains: 
2k B T . , ,,7. 

r=— smh^), ( 2 ) 

which in the limit of small (j/B) reduces to Newton's law for viscous deformation kinetics. 
We shall consider a steady transport of the glassy (or amorphous) electronic fluids in a wavy- 
rough microannulus of r\ (mean-averaged inner radius) with the inner interface being a fixed 
wavy-rough surface : r = n + e sin(/c# + (3) and r2 (mean-averaged outer radius) with the outer 
interface being a fixed wavy-rough surface : r = r<i + esin(/c#), where e is the amplitude of the 
(wavy) roughness, (3 is the phase shift between two boundaries, and the roughness wave number 
: k = 2tt/L (L is the wavelength of the surface modulation in transverse direction). 
Firstly, this amorphous fluid [6,12-13] can be expressed as 7 = 70 sinh(r/ro), where 7 is the shear 
rate, r is the shear stress, and 70 = B is a function of temperature with the dimension of the 
shear rate. In fact, the force balance gives the shear stress at a radius r as r = —[r 5(p e E z )]/2, 
where r is the shear stress along the boundaries of a control volume in the flow direction (the 
same direction as E z which (the only electric field) is presumed to be a constant or uniform), p e 
is the net charge density, | — S(p e E z )\ is the net electric force along the transport (or tube-axis 
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: z-axis) direction. 

Introducing \ = — ( r 2/2T~o)5(p e E z ) then we have 7 = 70 sinh(xr/r2). As 7 = —du/dr (u is the 
velocity of the (electronic) fluid transport in the longitudinal (z-)direction of the microannulus), 
after integration, we obtain 

u = u s H coshx — coshf — J , (3) 

X r 2 

here, u s is the velocity over the (inner or outer) surface of the microannulus, which is determined 
by the boundary condition. We noticed that a general boundary condition [13] was proposed 
for transport over an interface as 

An = L° 7 (l-^)- 1/2 , (4) 

7c 

where An is the velocity jump over the interface, L® is a constant slip length, 7 C is the critical 
shear rate at which the slip length diverges. The value of 7 C is a function of the corrugation of 
interfacial energy. 

With the boundary condition from [13], we can derive the velocity fields and electric-field-driven 
volume flow rates along the wavy-rough microannulus below using the verified boundary per- 
turbation technique [6,13-14]. The wavy boundaries are prescribed as r = r 2 + esin(/c#) and 
r = ri + esin(fe0 + /?) and the presumed steady transport is along the z-direction (microannulus- 
axis direction). 

Along the outer boundary (the same treatment below could also be applied to the inner bound- 
ary), we have 7 = (du)/(dn)\ on sur f ace - Here, n means the normal. Let u be expanded in e : 
u = uq + eu\ + e 2 U2 + ■ • •, and on the boundary, we expand u(ro + edr, 6(= 0q)) into 

u ( r ^)\(r +edr,9 ) = u(r Q ,0) + e[dr u r (r , 0)} + e 2 [—u rr (r o ,0)] H = 

7^2 X r 

{usiip + — [coshx - cosh(— )]}| surface , r = ri,r 2 ; (5) 

X r 2 

where 

u slip\on surface = on surface- (6) 

7c 

Now, on the outer interface (cf. [13-14]) 

du V(r — V2 — esin(/c#)) 

7 = Thl = |V(r - r 2 - esin(fc0))| ' ^ ' 

Considering ~ r\,r 2 S> e case, we also presume sinhx Tc/To- With equations (1) and 
(5), using the definition of 7, we can derive the velocity field (u) up to the second order : u(r, 0)= 
-(^27o/x){cosh(xr/r 2 )-coshx [l+e 2 X 2 sin 2 (/c0)/(2r|)]+exsinhx sm(k0)/r 2 }+u s ii p \ r=r2+esh ^ k e). 
The key point is to firstly obtain the slip velocity along the boundaries or surfaces. After lengthy 
mathematical manipulations, we obtain the velocity fields (up to the second order) and then 
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we can integrate them with respect to the cross-section to get the volume flow rate (Q, also 
up to the second order here) : Q = Jq p / r 7+/ B Tn(fce < + / 9) u ( r ' ®) r drdO = Q sUp + e Q Pl + e 2 Q P2 . In 
fact, the approximately (up to the second order) net electric- field-driven volume flow rate reads 
Q = Qout — Qin which is the flow within the outer (larger) wall : Q ou t without the contributions 
from the flow within the inner (smaller) wall Qi n : 

Q = 7r 7 o{L°(rf - rf) sinh X (1 - + ^[(4 - r\) cosh X - \r\ sinh X - 

7c/7o X X 

. . , ri. 2r 2 , , . , ri,„ 2r . , coshX/ r?., 

nr 2 smh(x— ) + -^-(coshx - cosh(x— )Jj + e {tt7 [x — ^ — (?"2 - — )J+ 

^ ■ • i /-, sinhx w i2 2 ,ru 2 . 7r r/ 
-7osmh X (l + --^)(-A; 2 + X 2 )[l - - 2 + ? 
4 7c/7o ^2 ^ 

^-^■coshx) + 7or2(-sinhx + X ) + 7o(n sinh(x — ) - — cosh(x— ))] + 
X X r 2 X r 2 

coshx, sinhx Nir 

7T7o{[sinhx + X (1 + . , . W2 - n cos/?J}+ 

fi 7c/7o 

TT 2 . COshx ri 2 

7 X 70— -j— [1 - (— ) ] cosh X . (8) 

4 7c/7o r 2 

Here, 

^ P() =L0 7 o[sinh X (l-|^)-V 2 ]. (9) 

3 Results and Discussions 

We firstly check the roughness effect upon the shearing characteristics because there are no 
available experimental data and numerical simulations for the same geometric configuration 
(microscopic annuli with wavy corrugations in transverse direction). With a series of forcings 
(due to imposed electric fields) : x = r 2[— S(p e E z )]/(2To), we can determine the enhanced shear 
rates (dj/dt) due to forcings. From equation (7), we have (up to the first order) 

'tL = ^°[sinhx + esin(^)^coshx]. (10) 
dt dt r 2 

Furthermore, if we select a (fixed) temperature, then from the expression of To, we can obtain 
the shear stress r corresponding to above forcings (x) : 

—i X 
T = Tosinh - [sinh(x) + e s'm(k9) — cosh(x) . (11) 

r 2 

Note that, based on the absolute-reaction-rate Eyring model (of stress-biased thermal activa- 
tion), structural rearrangement is associated with a single energy barrier (height) AE that is 
lowered or raised linearly by a (shear) yield stress r. If the transition rate is proportional to the 
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plastic (shear) strain rate (with a constant ratio : Co; 7 = CoRt, Rt is the transition rate in the 
direction aided by stress), we have 

where zv is an attempt frequency or transition rate, Co^o ~ 70 exp(Ai?//c^T), or 

T = 2 ¥Gk exp(A£/feT> if 5? <<L (13> 

The nonlinear character only manifests itself when the magnitude of the applied stress times the 
activation volume becomes comparable or greater in magnitude than the thermal vibrational 
energy. 

The following figures are for our primary interest : Higher (electric) conductivity or lower viscos- 
ity which means rather low resistance of the electric-field-driven transport of amorphous electron 
liquid. This can be manifested from equation (8) for the net transport Q (oc sinhx). We have 
the average velocity (for charge carriers) : v = Q/A m with A m being the effective area. After 
this, we then have the electric flux : p e v = |J| which conventionally equals to oE z where a is 
the electrical conductivity. Thus we can obtain the electrical resistance pn = 1/a. The latter is 
directly relevant to the frictional resistance or shear stress along the tube wall considering the 
balancing of the driven electrical force when the transport is finally steady or fully developed. 
This can be evidenced easily once \ IS small then we have p^ 1 oc v oc sinhx ~x a l _ 8(PeEz)\- 
Fig. 1 shows that there will be nearly frictionless or almost zero-resistance states if we select 
the activation energy to be 10~ 22 J ~ 0.001 eV. We can observe a sudden drop of the resistance 
(frictional or shear stress) around 3 orders of magnitude at T = 0.12°K (Vh ~ 10~ 24 m 3 ). It 
means there is a rather high electrical conductivity (or very low electrical resistivity) around 
this temperatrure for the material parameters selected. The qualitative as well as quantative 
similarity are that this critical temperature (T c phenomenon) resembles that found in amorphous 
superconductor Mg7oZn3o (cf. [9]). 

Similarly, Fig. 2 shows that there are almost zero-resistance states if we select the activation 
energy to be 10~ 22 J ~ 0.001 eV. We can observe a sudden drop of the resistance (frictional or 
shear stress) around 3 orders of magnitude at T = 1.9°K (Vh ~ 10~ 23 m 3 ). It means there is a 
rather high electrical conductivity (or very low electrical resistivity) around this temperatrure 
for the material parameters selected. The qualitative as well as quantative similarity are that 
this critical temperature (T c phenomenon) resembles that found in amorphous superconductor 
Fe x Nii_ x Zr 2 (x ~ 0.4; cf. Table 1 in [10]). 

To verify more of our approach, we demonstrate that if we select the activation energy to be 
10~ 22 J ~ 0.001 eV we can then observe a sudden drop of the resistance (frictional or shear 
stress) around 1 order of magnitude at T = 4.56°K (Vh ~ 10~ 23 m 3 ). It means there is a rather 
high electrical conductivity (or very low electrical resistivity) around this temperatrure for the 
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material parameters selected. The qualitative as well as qualitative similarity are that this criti- 
cal temperature (T c phenomenon) resembles that found in amorphous superconductor Zr7sRh25 
(cf. Fig. 1 (a) and Fig. 2 for the zero (magnetic) field case in [8]). 

Furthermore, Fig. 4 illustrates that there are almost zero-resistance states if we select the ac- 
tivation energy to be 3.5 x 10 -19 J. We can observe a starting (sudden) drop of the resistance 
(frictional or shear stress) at T ~ 164°K (Vh ~ 10 _20 m 3 ) around 3 orders of magnitude and be- 
low T = 80° K. This also demonstrate that there is a rather high electrical conductivity (or very 
low electrical resistivity) around this temperatrure for the material parameters selected. The 
qualitative as well as quantative similarity are that this critical temperature (T c phenomenon) re- 
sembles that found in amorphous high-T c superconductor Hgi_; C Ph r Ba2Ca2Cu308+ ( 5 [Hg(xPb) 
1:2:2:3] [15]. 

To give a brief summary, we already used the verified approach by borrowing the analogy from 
the successful application to the study of the almost frictionless transport of glassy solid he- 
lium to obtain the critical transport of electric-field-driven electronic fluids in metallic glasses 
as well as amorphous (bulk) superconductors. The critical temperatures we found after tuning 
of the material parameters resemble qualitatively as well as quantitatively those superconduct- 
ing temperatures reported before in (3D) amorphous superconductors. Our results show that 
effective shear-thinning is an important way to reach high-temperature charged superfluidity or 
superconductivity. We shall investigate other interesting issues [5,16] in the coming future. 
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Fig. 1. Calculated (shear) stresses or resistance using an activation energy 10~ 22 J or 
~ 0.001 eV. There is a sharp decrease of shear stress around T ~ 0.12°K. Around 0.12 K, 
the electric-field-driven transport of (glassy) electronic fluid is nearly frictionless. 
This critical temperature resembles that found in amorphous superconductor Mg7oZn3o 
(cf. [9]). 
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Fig. 2 Calculated (shear) stresses or resistance using an activation energy 10~ 22 J or 
~ 0.001 eV. There is a sharp decrease of shear stress around T ~ 1.9°K. Around 1.9 K, 
the electric-field-driven transport of (glassy) electronic fluid is nearly frictionless. 
This critical temperature resembles that found in amorphous superconductor Fe x Nii_ I .Zr2 
(cf. Table 1 in [10] : x ~ 0.4). 
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Fig. 3 Calculated (shear) stresses or resistance using an activation energy 10~ 22 J or 
~ 0.001 eV. There is a sharp decrease of shear stress around T ~ 4.56°K. Around 4.5 K, 
the electric-field-driven transport of (glassy) electronic fluid is nearly frictionless. 
This critical temperature resembles that found in amorphous superconductor Zr7sRh25 
(cf. Fig. 1 (a) and Fig. 2 in [8] : the zero (magnetic) field case). 



13 



0.7 
0.6 
0.5 



0.4 

03 
W 
CD 

3 0.2 



0.1 



Activation energy = 3.5 x 10" 19 Joule 



> 



Sudden drop 



Activation volume 



100 120 140 Q 160 

Temperature ( K) 



180 



Fig. 4 Calculated (shear) stresses or resistance using an activation energy 3.5 x 10~ 19 J. There is 
a sharp decrease of shear stress starting from around T ~ 164°K. Below around 80 K, the electric- 
field-driven transport of (glassy) electronic fluid is nearly frictionless. This critical temperature 
resembles that found in amorphous high-temperature superconductor Hgi_ a; Pb a ;Ba2Ca2Cu308_|_ ( 5 
[Hg(xPb) 1:2:2:3] [15]. 



